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■ Abstract. We study the ratio of bulk to shear viscosity in gluodynamics within a phenomenological quasiparticle model. We 

' show that at large temperatures this ratio exhibits a quadratic dependence on the conformality measure as known from weak 

coupling perturbative QCD. In the region of the deconfinement transition, however, this dependence becomes linear as known 
from specific strongly coupled theories. The onset of the strong coupling behavior is located near the maximum of the scaled 
1 interaction measure. This qualitative behavior of the viscosity ratio is rather insensitive to details of the equation of state. 
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INTRODUCTION 

' Transport coefficients, such as bulk (£) and shear (rj) viscosities, represent fundamental quantities specifying the 

physical properties of the matter. In the hydrodynamic description of the medium created in high-energy nuclear 
collisions they are important input parameters. First-principle lattice QCD calculations of the viscosity coefficients 
have been performed recently in SU C (3) gluodynamics in the deconfinement region [1, 2]. These lattice QCD results 
were shown to be successfully describable within a phenomenological quasiparticle model (QPM) for the gluon 
plasma [3, 4]. This model is based on an effective kinetic theory approach to QCD and includes non-perturbative 
effects via a thermal quasigluon mass. In particular, in the vicinity of the deconfinement transition temperature T c 
the shear viscosity to entropy density ratio exhibits a minimum [3] with a value as small as the predicted universal 
lower bound [5] rj/s ~ l/4n. This result is in clear contrast to naive extrapolations of standard perturbative QCD 
approaches into the non-perturbative regime close to T c . The specific bulk viscosity £/s is instead found to rapidly 
increase in the deconfinement transition region [3] . Moreover, the ratio 1^/77 is known to exhibit a distinct dependence 
on the conformality measure Av^ = 1/3 — with being the squared speed of sound. At weak coupling, i.e. at high 
temperatures T, the viscosity ratio depends quadratically on Av^ [6]. On the other hand, at strong coupling, i.e. close to 
T c , this dependence is expected to be linear [7]. The QPM reproduces the above asymptotic dependencies and exhibits 
a gradual interpolation between both regimes [4]. In this work, we analyze the structure of the viscosity ratio £/r] in 
SU C .(3) gluodynamics in the vicinity of T c , In particular, we discuss its sensitivity on details in the equation of state 
(EoS) and the corresponding v^. 
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THERMODYNAMICS AND TRANSPORT COEFFICIENTS 

In the QPM, the EoS, which describes changes of the pressure P with energy density e, as well as related ther- 
modynamic quantities are obtained from the local thermal equilibrium limit of the underlying effective kinetic the- 
ory. This model was shown to be very successful in describing the thermodynamic quantities obtained in lattice 
QCD [8]. Essential feature of the QPM is that the quasiparticle excitations follow a medium-modified dispersion 
relation, E 2 = p 2 +H(T). In a quasigluon plasma, the effective mass IT(r) = T 2 G 2 (T)/2 is quantified by the cou- 
pling G 2 (T) = 16n 2 / (ll\n{X(T-T s )/T c } 2 ), which near T c accommodates non-perturbative effects via only two 




FIGURE 1. Left panel: Scaled interaction measure (e — 3P)/T as a function of the scaled temperature T /T c for SU C (3) 
gluodynamics. The symbols depict lattice QCD results from [9] (squares), [10] (circles) and from [11] (diamonds). The dashed 
and solid curves are the QPM results for Fit 1 and Fit 2, respectively (see text). Right panel: Conformality measure Av^ as a 
function of T/T c in the QPM for Fit 1 (dashed curve) and Fit 2 (solid curve). 



parameters, A and T s . At large T, it reproduces the perturbative running coupling in QCD. To determine P and e in 
the QPM, one needs to introduce the integration constant B(T C ) as an additional parameter [8]. The squared speed of 
sound follows then from = dP/de. 

Figure 1 (left panel) shows comparisons of the QPM predictions for the scaled interaction measure (e — 3P)/T 4 as 
a function of T/T c with corresponding lattice QCD results [9, 10, 11]. The model parameters read as T s /T c — 0.73, 
A = 4.3 and B(T C ) = 0.19 T? (Fit 1) for data from [9, 10], while we use T s /T c = 0.52, A = 2.5 and B{T C ) = 0.48 T c 4 
(Fit 2) for data from [11]. The conformality measure obtained from these QPM results for the EoS is depicted in Fig. 1 
(right panel) as a function ofT/T c . The two different fits result in visible deviations in Av^ for T c <T < 2.5 T c . 

At the first-order phase transition at T = T c , the squared speed of sound v 2 (T) is discontinuous, which results in 
a discontinuity in Av 2 . We note that the limiting value A = \im T ^ T + Av 2 , which is linearly approached with T /T c as 
seen in Fig. 1 (right panel), is in general different from Av%(T c ) = 1/3. 

The bulk and shear viscosity coefficients follow directly from the effective kinetic theory [3, 12]. Their ratio was 
found in the following form [4] 
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where = T 2 ^frT 2 J^b / and = (r 2c j^j r 4 ^/^-2 depend non-trivially on T. The momentum integrals 

y k in £/ oa read as J k = J -0jn(T) [1 +d~ i n{T)] p 2 ~ k , where n(T) = d (exp(E/T) 1) 1 is the Bose distribution 
function for gluons with d = 16 and z denotes the relaxation time. For T > T c , these integrals follow the hierarchy 
^-2/T 2 J?o^> T 2 J^2 > 0, whereas dG 2 /dT 2 < 0. Assuming a momentum independent T, the ratio ^/r] in Eq. (1) 
is solely determined by parameters adjusted to equilibrium thermodynamics. 

In Fig. 2, the ratio C,/r\ from Eq. (1) is quantified for Fit 1 and 2. For both QPM fits one observes that the viscosity 

ratio for Av^ < 0.07, i.e. for T > l.5T c , is entirely determined by 15 (Av^J , i.e. by the quadratic, jz^-mclependent 
term in Eq. (1). This is a direct consequence of the behavior of the conformality measure at large T, which in the QPM 
reads as Av 2 ~ -75/ (18 dn 2 )T 2 (dG 2 /dT 2 ) + 0{G 2 T 2 (dG 2 /dT 2 )) with \T 2 (dG 2 /dT 2 )\ «G 2 «1. Thus, at leading 
order, all terms in Eq. (1) depend quadratically on Av^, where the first, ^ 2 -independent term dominates numerically. 

With increasing Av 2 , the £/rj ratio is reduced compared to 15 (Av;:) 2 by the non-perturbative terms in Eq. (1) which 
are proportional to (dG 2 /dT 2 ). Consequently, one finds ^/rj < 1 for T — > T c + . For Av^ > 0.17, i.e. for T < 1.15r c , a 
linear dependence on Av^ develops as seen in Fig. 2. In fact, near T c the factors l-sz^l become large as a consequence 

of a large \dG 2 /dT 2 \. This results in a cancelation of the quadratic dependence 15 (Av^) 2 in Eq. (1) by all other 
i*{),2-dependent terms. Effectively, Eq. ( 1 ) sums up near T c to a linear dependence of the form £/tj = a Av 2 + /3 . The 
transition to this linear behavior, as it is known for specific strongly coupled theories [7], takes place at T in the vicinity 
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FIGURE 2. Bulk to shear viscosity ratio £ /rj from Eq. (1) as a function of the conformality measure Av^ from Fig. 1 (right panel) 
for Fit 1 (long-dashed curve) and Fit 2 (solid curve). The short-dashed curve exhibits the quadratic dependence £/tj = 15 (Av,) 
and the dash-dotted curve shows a linear fit £/r/ = a Av^ + jS of the QPM result from Fit 1 with a — 3.78 and /3 = —0.305. For the 
QPM result from Fit 2, a similar linear fit yields a = 3.63 and /3 = —0.305. 



of the maximum in the scaled interaction measure, cf. Fig. 1. As is evident from Fig. 2, this qualitative behavior is a 
common feature of the viscosity ratio within the quasiparticle model for the gluon plasma irrespective of details in the 
EoS near T c . In the interval 0.07 < Av^ < 0.17, a gradual change between quadratic and linear dependence on Av^ is 
observed. We note that at T = T c the £/rj ratio develops a discontinuity in line with Av^. 



SUMMARY 

Within a phenomenological quasiparticle model for the gluon plasma, the bulk to shear viscosity ratio is found to 
exhibit, at large T, the quadratic dependence on the conformality measure as known from perturbative QCD. In the 
deconfinement transition region, this dependence becomes linear as found in specific strongly coupled theories. Thus, 
the quasiparticle model provides a systematic link between both regimes. The onset of the strong coupling behavior 
is located near the maximum in the scaled interaction measure. This qualitative behavior of the viscosity ratio is 
insensitive to details in the equation of state. 
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